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We examine hilltop quintessence models, in which the scalar field is rolling near a local maximum 
in the potential, and w « — 1. We first derive a general equation for the evolution of <f> in the limit 
where w «s — 1. We solve this equation for the case of hilltop quintessence to derive w as a function 
of the scale factor; these solutions depend on the curvature of the potential near its maximum. 
Our general result is in excellent agreement (Sw < 0.5%) with all of the particular cases examined. 
It works particularly well (8w < 0.1%) for the pseudo-Nambu-Goldstone Boson potential. Our 
expression for w(a) reduces to the previously-derived slow-roll result of Sen and Scherrer in the 
limit where the curvature goes to zero. Except for this limiting case, w(a) is poorly fit by linear 
evolution in a. 



I. INTRODUCTION 

Observational evidence p], suggests that approxi- 
mately 70% of the energy density in the universe is in the 
form of an exotic, negative-pressure component, dubbed 
dark energy. (For a recent review, see, e.g., Ref. [H). 
The observational bounds on the properties of the dark 
energy have continued to tighten. Taking w to be the 
ratio of pressure to density for the dark energy: 
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recent observational constraints are typically 
—1.1 w — 0.9 when w is assumed constant 
(see, e.g., @, [1] and references therein). 

One possibility is that the dark energy is, in fact, 
merely a cosmological constant, with w exactly equal to 
— 1. This is the standard ACDM model, in which the 
universe today is roughly 30% matter, and 70% vacuum 
energy, with the latter having a constant density. In this 
case, we expect future observations to continue to con- 
verge on this value of w. A second possibility is that w is 
close to —1, but not exactly equal to it. In this case, an 
alternative explanation for the dark energy is required. 

One possible model, dubbed quintessence, assumes 
that the dark energy arises from a minimally coupled 
scalar field; these models have been extensively explored. 
(See, for instance, Refs. [ESS IS EH for some of the ear- 
liest discussions). The equation of motion for the scalar 
field simplifies considerably in the limit where the back- 
ground expansion is well-described by ACDM, which is 
the case whenever w for the quintessence field is always 
close to —1. This fact was exploited in Ref. [ll[, which 
considered quintessence models in a potential satisfying 
the slow-roll conditions, i.e., 



approaches a single functional form, given by the present- 
day values of w and f2^. (Note that the term "slow- 
roll quintessence" refers to the fact that the potential 
satisfies the slow-roll conditions; the well-known slow-roll 
approximation for inflation cannot be applied [llf. For 
alternative approaches to this problem, see 

El)- 

However, the slow-roll conditions given by Eqs. ([2]) 
and ([3|), while sufficient to give w near —1 today, are not 
necessary. In this paper, we consider a second possibility: 
a scalar field rolling down near a local maximum in the 
potential. We assume that we are sufficiently close to 
the maximum that Eq. (|2|) still applies, but we relax our 
assumption in Eq. ([3]). In analogy with recent discussions 
of similar models in inflation [lfl EH E3 , we call these 
models "hilltop quintessence." A particularly important 
model that is well-described by this methodology is the 
pseudo-Nambu-Goldstone Boson (PNGB) model [II]. 

In the next section, we derive the version of the scalar 
field equation of motion that applies in the limit where 
the background expansion is well-approximated by a 
ACDM model. In Sec. Ill, we expand the potential 
for hilltop quintessence as a quadratic function near its 
maximum and solve this equation to derive the evolution 
of the scalar field and thus, the evolution of its equa- 
tion of state, for generic hilltop quintessence models. We 
find that all of these models can be characterized by a 
single set of functional forms for w(a) that depends on 
the present-day values of fl^ and w for the quintessence, 
along with a parameter (K) that depends on the curva- 
ture of the scalar field potential at its maximum. This 
result is given in Eq. (|3ip . Our results are discussed in 
Sec. IV. 
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It was shown in this paper that the evolution of w(a) in 
these models, which we will call "slow-roll quintessence," 



II. SCALAR FIELD EVOLUTION IN THE 
ACDM LIMIT 

Consider a minimally-coupled scalar field, </>, in a po- 
tential V((f>). The density and pressure of the scalar field 
are given by 



V{<j>), 
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and 
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respectively, and the equation of state parameter, w, is 
given by equation |T]). The equation of motion for this 
field in an expanding background is 



• dV 
■ + 3H(/>+ — = 0, 



(6) 



We expect Eq. (jl3|) to provide a good approximation 
to the evolution of any scalar field in the limit where 
the Hubble parameter in Eq. (J6j) is well-approximated 
by ACDM, i.e., any model in which the kinetic term in 
Eqs. (|4]) and ([5]) is dominated by the potential term, so 
that w for the scalar field never evolves very far from 
— 1. Eq. (I13|) applies, for example, to the limiting model 
examined in Ref. [TlJ. We now apply it to a different 
class of models. 



where the Hubble parameter H is given by 

h=(-) = V^A 



(7) 



and we have assumed a flat universe. Here a is the scale 
factor, pt is the total density, and we work in units for 
which 8ttG = 1. The evolution of the scale factor is 
further described by 
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where pr is the total pressure. 

We first eliminate the <fi term in Eq. 
change of variables 



4>{t) = u{t)/a{tf/\ 
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by making the 
(9) 



which gives 
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u + a 3 ' 2 V l {u/a 3 ' 2 ) =0. (10) 



Applying the expressions for a/a and a/a from Eqs. 
and © gives 



u+ -p T u + a 3/2 V'(u/a 3/2 ). 



(11) 



This equation is as yet exact. 

Eq. (|11[) takes a particularly simple form for two 
special cases: a matter-dominated universe (for which 
Pt = 0) and a ACDM universe (for which pt is con- 
stant). In the former case, we obtain 



u + tV'(u/t) = 0, 



(12) 



and it is easy to use this equation to derive, in a straight- 
forward way, the results of Refs. [6] and [9-]. However, 
the subject of this paper is the second case. 

We assume a universe containing matter and 
quintessence, but such that the quintessence always has w 
near —1. In this limit, pr is well-approximated by a con- 
stant: pt ~ —p<f>Q, where p^ is the nearly constant den- 
sity contributed by the quintessence in this limit. Then 
Eq. (fTTj) becomes 



\p^u + a 3 / 2 V'{u/a 3 ' 2 ) = 0. 



(13) 



III. HILLTOP SCALAR FIELD EVOLUTION 

We are interested in models in which the scalar field 
evolves near a local maximum in the potential. The most 
important model of this type is the PNGB model [18j , for 
which the potential is given by 



V(cf,) =Af 4 [cos(0//) + l], 



(14) 



where M and / are constants. (For recent discussions 
of the PNGB model in the context of dark energy, see, 
e.g., Refs. [H, H3, HII and references therein). Other, 
less- well-motivated models with a local maximum in the 
potential include the Gaussian potential, 



V(4>) = M 4 
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and the quadratic potential 

v(4>) = v -v 2( j ) 2 . 
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Our purpose in this paper, however, is not to examine an 
exhaustive list of such models, but to show that they all 
converge to a common evolution under certain conditions. 

We note that any model with a local maximum in the 
potential at <j) = 0* can be expanded, near this maxi- 
mum, in the form 



V(<f>) = V(<t>*) + (l/2)y"(0*)0 2 + O(0 3 



(17) 



This expansion will be a good approximation for the 
PNGB model when <j> ^ f and for the Gaussian model 
as long as 4> o. It is exact for the quadratic potential. 
The evolution of <f> for this expansion (or, equivalently, 
for the quadratic potential of Eq. I16p was previously 
examined in Refs. [13, [H, H3, H3l- The motivation in 
Ref. 1221 was similar to our own investigation, while 
Refs. |23l I24I [25l ] were concerned with the future fate 
of the universe in such a model. The exact solution to 
the scalar field equation of motion for this potential is 
given in Ref. 12211 for the matter-dominated case, while 
Refs. [H, [24], |25| | give the evolution for the scalar-field- 
dominated case (appropriate for the far future of the uni- 
verse). Here we solve for the general case in which both 
matter and the scalar field contribute to the evolution, 
since this is the relevant case at low redshift. 

Substituting the expansion of Eq. (fT?]) into Eq. (TTT 
(and taking = ^ (</>*)) gives 



u + [V"(4>*) - (3/4)V(<M]« = 0. 



(18) 
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If we define the constant k to be given by 



k = V(3/4M0*) ~ V" {</>.), (19) 

then the general solution to Eq. (fT8| is 

u = Asinh(fci) + Bcosh(fci). (20) 

Now we assume that the scale factor is well-approximated 
by its value in the ACDM model, which is again a con- 
sequence of the assumption that w is always close to —1: 
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where Sl^o is the present-day value of and a = 1 at 
present. The time t\ is defined to be 



(22) 



With this expression for a(t), the general solution for (f){t) 
is given by 



n 



-,1/2 



A sinh(fct) + B cosh(fci) 
sinh(f/*A) 



We require that at t = 0, is equal to a fixed initial 
value, 0j. This forces B = and gives us the value of A, 
so that 



sinh(fci) 



kt A sinh(i/tA) 



(24) 



In the limit where t <C t A , our solution reduces to the 
matter-dominated solution in Ref. [22j|, while for t 3> 
iA, we regain the scalar-field-dominated solution of Refs. 

The equation of state parameter for quintessence is 
given by 



1 + w 



(25) 



(23) Taking p « p 0o w V(<£„), Eqs. (|24j) and (25]) yield 
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We can normalize this expression to the present-day value express wasa function of the scale factor (or the redshift) 
of w, which we denote wq, and we can use Eq. (|2ip to rather than t. We obtain: 



1 + w(a) = (1 + w )a 



'TT^ ) kt A cosh [kt (a)] — — Clfoja 3 + f20osinh [kt (a)] 



[v/O^o^A cosh(fcio) — sinh(fcio) 
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where t(a) and to can be derived from Eq. (j2Tj) : 
t(a) = ^a sinh -1 

and 
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i = tAtanh 1 (Vf^o) 



(28) 
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We now define the constant X = ht A . In terms of the 
quintessence potential, K is just 



K = y/\-{A/Z)V"{^)/V{<t>,). (30) 
Thus, K depends only on the value of the potential 
and its second derivative at its maximum. Note that 
V" ((/>*) < 0, so K > 1. Then Eq. (gTJ) can be written as 



+ »(„) = (1 + + l)K{K - F [f + iFi t 1)K<K + Fi ^S 



[(fi- 1/2 + l)^(K-O7 1/2 ) + (f77 1/J 



(31) 
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FIG. 1: The evolution of w(a) given by Eq. (pi) for hilltop 
quintessence models with wo = —0.9 and Q^o = 0.7, for the 
indicated values of K, as defined in Eq. (|30[) . Red curve (filled 
circles) gives the approximation from Ref. for slow-roll 
quintessence. 



where F(a) is given by 
F(a) = 



1 + 



Da- 



rn 



(Note that F(a) = 1/ \/fl^(a), where f2^(a) is the value of 

) 



1) 
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fl^ as a function of redshift, so that F(a 

Eq. (|3ip is our main result. It shows that, in the limit 
where w is close to —1 (i.e., the scalar field potential 
energy dominates the kinetic energy), all of the hilltop 
quintessence models with a given value of wq form a sin- 
gle set of models parametrized by the value of K, which 
depends only on V"(fa)/V(fa). For the case of slow-roll 



quintessence models, a similar analysis yields only a sin- 
gle form for w(a), once Wq and fl^o are fixed [11] ■ Here 
we have more complex behavior, since w(a) also varies 
with K. This behavior is illustrated in Fig. 1, where we 
fix Wq — —0.9 and f^o = 0.7 and then plot w(a) from Eq. 
(f3TJ) as a function of a. When V" (4>o) /V (4> ) > 1, we see 
that the field evolution diverges significantly from slow- 
roll quintessence. In particular, w increases more slowly 
at high redshift, but much more rapidly at low redshift. 
This is easily understood from the nature of these po- 
tentials: dV/d(f> increases as <fi rolls down the potential. 
In the limit where V" (<po)/V((j)o) — > 0, corresponding to 
K — ► 1, both slow-roll conditions (Eqs. [2] and [3]) ap- 
ply, and we expect the evolution to converge to the form 
derived in Ref. [11], namely 



1 + W = (1+UJo)' 
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(33) 

It is straightforward to show that Eq. (f3"Tj) reduces to 
Eq. (|33"|) in the limit where K — > 1. We also show this 
effect in Fig. 1: the curve corresponding to K — 1.01 is 
indistinguishable from w(a) for slow-roll quintessence. 

Eq. (f3"Tj) simplifies considerably for the case where K 
is a small integer. We have, for example, 
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Of course, it would require rather bizarre fine-tuning of 
the potential for K to be exactly equal to one of these 
values, but these special cases provide some qualitative 
insight into the behavior of w(a) as a function of K. 



Now we evaluate the accuracy of Eq. (|3T|) when applied 
to various models of interest. Consider the three hilltop 
potentials outline d above. We hav e, for the quadratic 
potential, K = yjl + (8/3)(V 2 /V Q ), for the PNGB po- 
tential, K = (2/3) (l// 2 ), and for the Gaussian 
potential, K = y/l + (8/3)(l/cr 2 ). In Figs. 2-4, we have 
plotted our expression for w(a) from Eq. (pTTj) against 
the exact numerical evolution for these three potentials, 
fixing Qj,q — 0.7. (Once the potential and fJ^o are fixed, 
the value of fa is chosen to produce the desired value of 

W ). 



The agreement, in all three cases, between Eq. (|3"Tj) 
and the exact numerical evolution is excellent, with er- 
rors Sw < 0.5% for the quadratic potential, Sw < 0.1% 
for the PNGB potential, and Sw < 0.3% for the Gaus- 
sian potential. Rather surprisingly, the errors for the 
quadratic potential (which is the basis of our approxima- 
tion) are actually larger than for the other two potentials. 
However, this is due to the fact that we have made several 
approximations in our calculation; these tend to cancel 
for the PNGB and Gaussian potentials. 

The agreement between Eq. (|3"Tj) and the exact evo- 
lution is particularly striking for the case of the PNGB 
potential, which also happens to be the most interesting 
and well-motivated of the hilltop potentials. We have 
therefore examined this potential in more detail. In Fig. 
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FIG. 2: Comparison between our approximation for w(a) for 
hilltop quintessence (Eq. 13 1 p with wq — —0.9 and fi^o = 0.7 
and the exact (numerically-integrated) evolution for w(a) for 
the quadratic potential V — Vo — Vify 2 . Red (dashed) curves 
give our approximation, and solid (blue) curves give exact 
evolution, for (left to right), K = 1.01,2,3,4, where K is 
defined by Eq. ([50)1. 
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FIG. 3: As Fig. 2, for the PNGB potential, V = 
M 4 [cosW>//) + l]. 



5, we extend our results to larger values of wq. It is in- 
teresting to see that the agreement remains excellent for 
wq as large as —0.7. Thus, Eq. (|3"Tj) represents a nearly 
exact solution for w(a) for the PNGB model within a 
wide range for wq. 

Finally, in Figs. 6-9, we use Eq. (|3ip to construct a x 2 
likelihood plot for wo and f^o with K = 1.01,2,3,4, 
using the recent Type la Supernovae standard candle 
data (ESSENCE+SNLS+HST from [5]). While none of 
these models is ruled out by current supernova data, it 
is interesting to note that the hilltop quintessence mod- 
els (K = 2, 3, 4) produce a larger allowed region than the 
slow-roll quintessence model (K — 1.01), and the allowed 
region increases with increasing K. 



FIG. 4: As Fig. 2, for the Gaussian potential, V = 

M 4 e-^ 2 . 
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FIG. 5: As Fig. 3, for the PNGB potential, extended to larger 
values of wo. 



IV. DISCUSSION 

Our results indicate that hilltop quintessence models 
with w near —1 all produce a similar evolution for w(a), 
given by Eq. (|31[) . The importance of this result lies in 
the fact that a very general set of models can be mapped 
onto a fairly constrained set of behaviors for w(a). Note 
that in general, the evolution given by Eq. (|3 1 1) and 
shown in Fig. 1 is not well-described by the popular 
linear parametrization, w(a) = Wq + w a (l — a) [26j |. The 
one exception is the limiting case K — * 1, where we regain 
the results of Ref. |ll|: in this limit the evolution is 
roughly linear for a > 0.5 {z < 1), with w a ~ —1.5(1 + 
w ). 

This investigation, along with Ref. [ll|, can also be 
thought of as a kind of Taylor expansion of the poten- 
tial about the initial value of <f>. The results of Ref. [ll| 
apply when the linear term dominates, while the results 
presented here assume a quadratic expansion. The re- 
sults of Ref. [TTf| are then a special case (albeit a very 
important special case) of the results presented here. 

Of course, our results do not apply to all thawing 
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FIG. 6: Likelihood plot from SNIa data for the parameters wo 
and Sl^o, for hilltop quintessence models with generic behavior 
described by Eq. (|31[) . with K = 1.01, where K is the function 
of the curvature of the potential at its maximum given in Eq. 
(|30|) . The yellow (light) region is excluded at the 2a level, and 
the orange (darker) region is excluded at the la level. Red 
(darkest) region is not excluded at either confidence level. 



FIG. 9: As Fig. 6, for K = 4. 




quintessence models with w near —1, but only those sat- 
isfying Eq. ([2]). A different approach was taken recently 
by Cahn, et al. [3| • They looked at the evolution of the 
scalar field at early times, when pt S> p<f>, while assuming 
nothing about the detailed nature of the potential. This 
allows the derivation of a generic result for the evolution 
of the scalar field before it dominates the expansion (see 
also Ref. [27| for the special case of tracker fields). Our 
results for the particular class of potentials considered 
here agree with those given in Ref. [l4| (as they should) 
in the limit where Qh, <C 1. 
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FIG. 7: As Fig. 6, for K = 2. 




FIG. 8: As Fig. 6, for K = 3. 
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